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We analyze the instanton transitions in the framework of the gauge invariant variational calcula¬ 
tion in the pure Yang-Mills theory. Instantons are identified with the saddle points in the integration 
over the gauge group which projects the Gaussian wave functional onto the gauge invariant physical 
Hilbert space. We show that the dynamical mass present in the best variational state provides an 
infrared cutoff for the instanton sizes. The instantons of the size p < 1/M are suppressed and the 
large size instanton problem arising in the standard WKB calculation is completely avoided in the 
present variational framework. 
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I. INTRODUCTION 

In recent years there have been renewed interest in application of the Hamiltonian methods to the study of non- 
abelian Yang-Mills theories ||l|-^ . One set of these works attempts to solve “exactly” the strongly interacting 

gauge theories in the sense that a nonlinear transformation is performed to a set of gauge invariant coordinates. 
One then tries to find a controlled expansion akin to strong coupling perturbation theory, which hopefully solves 
the infrared part of the theory in the leading order. Another set attempts to find the Yang-Mills vacuum wave 
functional with the help of a variational approximation. In particular in Q a gauge invariant generalization of a 
Gaussian variational approximation was developed. The hope of this approach is that the vacuum of QCD may be 
not very different from the vacuum of a free theory in many important respects. This hope rests on the observation 
that many genuine nonperturbative effects in QCD appear already on the momentum scales much larger than Aqcd 
where the coupling constant is still small j^. It may be possible then to account for these effects with the Gaussian 
wave functional, which is similar to the ground state of the free theory modifying only the width of the Gaussian for 
the low momentum modes. This modification is essentially nonperturbative and should lead to the generation of the 
same condensates that account for a variety of QCD phenomenology via the QCD sum rules [^. 

The variational calculations of are of exploratory nature and many questions regarding their validity remain 
unsettled. Some of those are discussed in the original paper Q as well as in ||^. Nevertheless this variational approach 
captures many of the essential features of gluodynamics; mass generation, formation of the gluon condensate, ||], and 
asymptotic freedom, ,|l^ . Possible appearance of the linear potential between static quarks in this approach 
has also been discussed 1 9[-|l 1|| . In fact one of the nicer features of this approximation is that it exhibits nontrivial 
nonperturbative infrared physics (gluon condensate) along with correct weak coupling ultraviolet behavior (one loop 
Yang-Mills /3 function). If so one is naturally lead to ask whether it also gives a good account of the instanton physics. 
Instantons are the only concrete description of the non-perturbative nature of the QCD vacuum in the path integral 
formalism. In the ultraviolet region although the effect of the instantons is nonperturbatively small they are easily 
identifiable. They should therefore serve as a useful probe for any non-perturbative approach especially if it purports 
to capture both the infrared and the ultraviolet physics. The aim of the present paper is precisely to study the 
structure and the properties of the instantons in the variational approach of ^ . 

Instantons are localized, finite-action classical solutions of the field equations of QCD in Euclidean space-time [ p^ . 
Such solutions have been obtained in exact analytical form, and extensively studied. An excellent review of instantons 
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in gauge theories can be found in |l^ and 0 - 

Physically instantons represent tunneling processes between topologically distinctive vacuum sectors with the ex¬ 
ponent of the instanton action being equal to the transition probability between two of these vacuum states. As any 
tunneling probability at weak coupling it is nonperturbatively small (of order exp(—const/as)) and the instantons 
therefore are invisible in the weak coupling perturbation theory. When first discovered there was hope that instantons 
would provide the solution to the strong coupling problem Although it has been subsequently realized that 

instantons are irrelevant for understanding confinement, they have provided a beautiful mechanism of spontaneous 
chiral symmetry breaking |p^ . The instanton liquid model (initially introduced on phenomenological grounds [p^ , 
and later justified by Euclidean variational methods |[l9|] ) to this day remains the most complete theory of chiral 
symmetry breaking in QCD. It is therefore vital to understand the properties of the instantons if one is hoping to 
extend the application of the Gaussian variational approximation to QCD with fermions. 

Although the notion of the instanton is intrinsically Euclidean, the tunneling between different vacuum sectors can 
be formulated in the Hamiltonian language as well as in the Lagrangian one. In fact the gauge projected formalism 
of 1^ is very well suited for this purpose. The projection of the initial Gaussian onto the gauge invariant subspace is 
achieved by the integration over the gauge group. As will be explained in detail later, the SU{N) matrix U{xi) (the 
dependence here is on the spatial coordinates only) in this type of calculation turns out to play the role of a cr-model 
field. This field is governed by an “action”, whose structure depends on the parameters of the variational state. As 
we shall see, the Yang-Mills instantons correspond to the topologically nontrivial saddle points of this action. In this 
paper we will study the properties of these saddle point solutions. 

Our main result is the following. We find that the appearance of the dynamical mass parameter in the variational 
state stabilizes the size of the instantons. Recall, that the main unsolved problem of the dilute instanton gas ap¬ 
proximation is that (when account is taken of the one loop running of the coupling constant) the path integral is 
dominated by the large size instantons. The measure for the integration over the sizes diverges as a power in the 
infrared and this divergence renders the dilute gas calculation meaningless. It is usually assumed that this infrared 
divergence is eliminated by some nonperturbative effects. This is precisely what we find in the variational approach. 
The dynamical mass M present in the best variational state suppresses instantons of sizes p > 1/M. The size of the 
stable instanton that we find turns out to be consistent with the average size of the instanton in the instanton liquid 
model. 

This paper is structured as follows. In section || we briefly recall the formalism and the results of . We explain how 
to identify the tunneling transition in this framework and what type of classical configuration should be identified 
with the instanton. It is also noted that the generation of the dynamical mass itself found in |^] can be directly 
interpreted in terms of the condensation of these instantons. 

In section III we study numerically the structure and the action of the small size instantons. For these instantons 
the presence of the dynamical mass is irrelevant and this calculation is performed at zero mass. The profile function of 
the instanton and the value of the transition probability are approximately determined by a variational method. We 
find that the tunneling probability indeed scales as exp{—ccts} with the value of c approximately two times larger than 
in the standard Euclidean calculation. We explain why this discrepancy is not unexpected in a variational calculation. 

In section IV corrections to the solution and the action due to a non-zero mass gap are calculated. It is found that 
the instanton is stabilized to a size of the order of the inverse mass scale. This size is directly comparable with that 
found in the instanton liquid model. 

Finally, section^ is devoted to discussion of our results and their relation with the instanton liquid model. 


II. THE HAMILTONIAN PORTRAIT OF AN INSTANTON 


We start with a brief description of the gauge invariant Gaussian approximation of Q. The Ansatz for the QGD 
wavefunctional considered in B is. 


with, 


vE-[A“] = J DU{x,)'l^u[A], 

^u[A] = exp |-i J SxSy A^%x)G-l<^\x - y) A; 


{y) 


Here A^ is the gauge transform of the vector potential with the gauge transformation U. 

At{x) = S‘^\x)Al{x) + Xt{x), 


( 1 ) 


( 2 ) 


(3) 
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with. 


S‘^\x) = itr [r^U^xyUix)] , (4) 

Xi{x) = -tr \T°'U^{x)diU{x)] . 

9 

The SU(N) generators are taken to satisfy the following algebra and normalization, 

itr[r“T'’] = (5) 

The state thus constructed obeys Gauss’ law since it is explicitly invariant under the gauge transformation Ai{x) 

Ay (x) with arbitrary SU{N) matrix V{x). The width of the Gaussian G{x) is a parameter with respect to which 
the expectation value of the Hamiltonian is varied. A simple rotation and (global) color invariant form is = 

. The functional form of G~^ is chosen to agree with the perturbation theory in the limit of high momentum 
on one hand, and to allow for the nonperturbative mass scale on the other, 


G-\k) 


if 

M if P < 


( 6 ) 


Technically, the calculation of the expectation values of gluonic operators in the state eq.(||) is mapped into the 
calculation in a nonlocal nonlinear tr-model in three Euclidean dimensions. Gonsider the vacuum average of an 
arbitrary gauge invariant operator 0[A], 


(O) = i J DA'f’*[A]0[AmA] 



j DUDA-^1[A]0 [A]^u[A]. 


DUxDU2DA-^*jj^ [A]0[A]-^U2 [^] 


(7) 


In the last line the matrix U = U\U 2 is the relative gauge transformation between the two Gaussian wave functions. 
The normalization factor Z (the norm of the state) is, 


Z = j DUDA-^\[A]^u[A]. (8) 

The Gaussian integration over the gauge potential Ai can be performed explicitly. As a result the last step of the 
calculation is a path integral over the SU(N) matrix U with the cr-model partition function, 


Z = J DU exp{-T[U]}, (9) 

where the (nonlocal) action is, 

r[17] = ^AAA + ^TrlnM. (10) 

Here summation over all indices (rotational, color and coordinate) is implied. The first term is written explicitly as, 

J Px Py \1{x)/S°‘^[x,y)Xl{y), (11) 


with. 


In eq. & we have defined, 


A“^(x,y) = [<5“'’G(x -y) + S<^\x)G{x - y)S<^\y)]-^. 


M^Hx,y) = + H G-\x - y) 5,,. 


( 12 ) 
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The second term in eq. (0) is of 0{g^) relative to the first one and with the accuracy of can be ignored. We will 
not consider it in the following. The action of the tr-model eq. (0) depends on the variational parameter M through 
eq. ©• The minimization of the expectation value of the energy in leads to a nonzero value of the mass parameter 
M in the best variational state. The dynamical mass parameter is determined by the relation 

as{Mo) = (13) 

where the Yang-Mills coupling constant Us evolves according to one loop /l-functionQ The significance of this value of 
M from the point of view of the effective a model is that at this point it undergoes the phase transition. For M < Mq 
the model is in the weakly coupled ordered phase. In this phase the matrix U is close to the unit matrix with small 
fluctuations around it. For M > Mq the model is in the disordered phase. The matrix U fluctuates strongly so that 
it covers all available phase space and its average value vanishes. It was found in Q that the energy is minimized just 
above the phase transition M = Mq+ so that the cr-model is in the disordered phase. 

Coming back to the subject of the present paper, the first thing is to understand how do we expect to see instantons in 
this formalism. The answer to this is the following. As explained above the effective cr-model arises as an integration 
over the relative gauge transformations between the two Gaussian states in the linear superposition eq.(|^). The 
Boltzmann factor exp(—r[17]) for a given matrix U is therefore just the overlap of the initial and the gauge rotated 
state, or in other words the transition amplitude between the two states^. The instanton transition is precisely the 
transition of this type, where the two states are related by a large gauge transformation. The matrix of this large 
gauge transformation must carry a nonzero topological charge Il 3 {SU{N)). 

The integration measure over U indeed includes integration over topologically nontrivial configurations. The finite¬ 
ness of the action eq. © requires that the matrix U approaches constant value at infinity. This identifies all points at 
spatial infinity hence, the physical space of the model is S^. Field configurations are maps from into the manifold 
of SU{N) and are classified by their winding number, or topological charge, which is an element of the homotopy 
group Il 3 {SU(N)) = Z. The tr-model action in a given topological sector is minimized on some configuration which 
is a solution of classical u-model equations of motion. In particular, the solution with a unit topological charge is 
expected to have a “hedgehog” structure much like the topological soliton in the Skyrme model j^. The integral 
over U in the steepest descent approximation is saturated by these classical solutions. 

These u-model configurations that belong to a nontrivial topological sector with a unit winding number represent 
QCD transitions between the topologically distinct sectors. The topologically nontrivial classical soliton solutions of 
the cr-model are therefore the three dimensional images of the QCD instantons. 

The QCD instantons are defined in space time and are therefore four dimensional point-like objects. The tr-model 
solutions are intrinsically three dimensional. Nevertheless, there is a natural simple relation between the two. For a 
given Yang-Mills instanton solution one can find a three dimensional SU{N) matrix U{xi) by the procedure 

discussed by Atiyah and Manton, pl| , 

UAM{xi) = Pexp , (14) 

where the contour of integration C is a straight line Xi = const^ —oo<Xq<oo. The matrix Uam gives the relative 
gauge transformation between the initial trivial vacuum at a:o ^ — oo and the topologically nontrivial vacuum at 
xq -foo, or in other words between the initial and final states of the instanton transition. Clearly, its meaning is 
precisely the same as of the classical soliton solution of the effective cr-model eq. ( 0 - Also, the QCD instanton action 
and the tr-model soliton action have the same physical meaning. They both give the transition probability between 
different topological sectors in QCD. We will therefore refer to the cr-model solitons as instantons in the following. 

One has to realize that although the QCD and the a - model instantons have the same physical meaning, it is 
not assured that the numerical value for their respective actions is the same. They both approximate the value of 
the transition probability in QCD, but the approximations involved are quite different. The QCD instanton action 
is the result of the standard WKB approximation which is valid at weak coupling and therefore for small instantons, 
but breaks down for instantons of large size. The cr - model instanton action on the other hand is the value of this 
transition probability in a particular Gaussian variational approximation. It is natural to expect that variational 


^More accurately, the /3-function in the variational calculation of Q is slightly different from the complete one loop expression. 
See discussion in section 0 - 

^Since the space of the matrices U is continuous, strictly speaking the Boltzmann factor is the differential rather than the 
total amplitude. 
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calculation underestimates the value of the transition probability at very weak coupling. The transition probability 
is given by the overlap of the “ground state” wave functions in two topological sectors. For simplicity let us consider 
a quantum mechanical system with two vacua at x±. If the area below the barrier separating the vacua is large, 
the standard WKB instanton calculation is applicable. The wave function of each of the vacua below the barrier 
has essentially an exponential fall off exp{i ^E — V{x — x±)}. The instanton calculation is the calculation of the 
overlap of these functions. Our variational calculation corresponds to approximating the respective “ground states” 
at x± by Gaussian wave functions. The tails of the Gaussians fall off much faster away from the minimum than the 
actual wave function and the overlap is therefore is expected to be smaller. When the coupling constant is not too 
small (or when the area below the barrier is not too large) the overlap between the two states is not determined any 
more by the behavior of the “tails” of the wave functions. In this situation one can expect the Gaussian approximation 
to do much better, since the overlap region contributes significantly to the energy and therefore plays important role 
in the minimization procedure. 

The rest of this paper is devoted to a quantitative study of instantons in the variational state eq.(|^). Before 
proceeding to this part, however, we would like to note that implicitly the instantons played a very important role 
already in the energy minimization of |^ . As we mentioned above the energy is minimized for the value of the mass 
parameter M at which the a model is in the disordered phase. The transition between ordered and disordered phases 
in a statistical mechanical system can usually be described as a condensation of topological defects. This is a standard 
description of the phase transition in the Ising and XY models [Q. In the a model eq. (0) the relevant topological 
defects are none other than the instantons. In this sense the appearance of the dynamical mass in the best variational 
state itself is driven by the condensation of instantons. 


III. SMALL SIZE INSTANTONS 

In this section we wish to study small size instantons. For the instanton solutions of a size p « 1/M the presence 
of a finite dynamical mass is irrelevant. We will therefore take M — 0 for the calculations in this section. The 
existence of a finite mass scale is very important for the instantons of large size and will be taken into account in the 
next section. It somewhat increases the complexity of the calculation but all the necessary methods can be developed 
for the M = 0 case. 

For M = 0 the inverse propagator (||) which defines the variational state eqs.(P|) is (in momentum space) G~^{k) = 
|fc|. We find it more convenient to work in coordinate space throughout the rest of this paper. Fourier transforming 
G~^ to the coordinate space according to the definition, 

G-\k) = (27r)-3/2 J d^xG-\x)E'^^, (15) 


we find. 


G ^{x-y) 


]_ ( 0(|a:-?/| - A ^) 
V \x-y\^ 


A^5(|a; 


2/1-A-i)). 


(16) 


Here we had to introduce the ultraviolet cutoff A to define the coordinate space expression properly. The coefficient 
of the second term is determined by requiring that at finite cutoff A, 


J dx^G ^{x)(xG ^(fc = 0) = 0. 


(17) 


Note that A is introduced as a regulator only and no subtraction in the expression for the propagator (or the action) 
was performed. The cutoff A should be taken to infinity at the end of the calculation and the results of the calculation 
should be finite in this limit. Below we show how the divergent terms cancel exactly, so that in the numerical 
calculations we only take into account finite terms in the limit A —> oo. 

We are searching for an instanton solution of topological charge one to the action &■ Such a solution should have 
the maximally symmetric “hedgehog” form, 

U{x) = e^'^ ^ = cos/(r) + ir“a;“sin/(r), (18) 


where r“ are the generators of SU{2) and /(r) is an unspecified function of r = |a;|, which we will refer to as the 
profile function. The profile function is constrained to satisfy /(O) = tt and /(oo) = 0, which ensures that the field 
configuration eq.(^ has unit topological charge. We need only consider the group SU(2) since the solutions for 
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SU{N > 2) can be found from the embedding of SU(2) in SU(N), ||^. Substituting this form for the held into the 
components of the action, we hnd, 


A“(a;) = {x)diU {x)] (19) 

9 

= --x^x°-f{r) - —((5“ - i*x“)sin2/(r) - -eabi—sin^ f{r), 
g gr 9 r 

S'^^ix) = ^tr[T‘^U+ixyU{x)] (20) 

= (5“** - eahci"" sin 2/(r) + 2(i“x'' - (5“^) sin^ /(r), 

where Cabc is the antisymmetric tensor. For the hedgehog conhguration, the function S°'^{x) differs signihcantly from 
6°"’’ only over a small region, the size of which depends on how fast the transition is from the asymptotic behavior at 
large distance to that at small distance. If the prohle function has a sharp transition between these two limits one 
can reasonably approximate 3°^^ by (5“**. Once the form of the prohle function has been found this can be checked 
for consistency. It is possible to expand around 3°^^ = 5“^ and systematically calculate corrections in 3°“^ — 6°“^. At 
the end of this section we calculate the leading correction and hnd that it is indeed small for the instanton prohle 
function. 

This approximation allows us to write. 


A“^(x,y) = [S-^Gix -y) + 3-\x)Gix - y)3-\y)]-^ ^ 
With the dehnition eq. (0), we write the action eq. (0) as, 


r = ri-r2. 


Ti 

r2 


1 

dTT^ 

47r^ 


d^xcPy\^{x) 


Qd^-yl 

\x-y\‘^ 




d^xSyXl{x)5{\x-y\-K ^)A“(y). 


( 21 ) 


( 22 ) 


Both Fi and r 2 are divergent. The divergences between the two terms however cancel leaving the action hnite. 
Changing variables, Ui = Xi + yi, Vi = Xi — yi, and using the fact that = V^|u|“^, we hnd after integrating by 

parts twice. 


Fi 


Tsurf 


Tdiv 


[A“(a:)A“(y)] + F.^rf + Fdiv, 
(2kr^ +5|^|) [A“(a:)A“(y)] ^ ^ 
f d^udriy {2\v\-^ + d\^i) [X'}ix)X'}{y)] 


(23) 

(24) 

(25) 


where d^v denotes angular integration in u-space, and Vj is the ?;-space Laplacian. The divergent part Fdiv cancels 
exactly the “subtraction” term r 2 when A —> oo. Moreover, the surface terms Fgurf vanish if the prohle /(r) approaches 
zero as or faster at inhnity, and is inhnite otherwise. In what follows we assume that /(r) decreases fast enough. 
Ignoring all terms in the action which vanish in the limit A ^ oo we are left with the hnite, cutoff independent action, 


1 


1 






52 


dxj dyj 


[K{^)K{y)]- 


For the hedgehog conhguration 


K{x)Xiiy) 


4cos2 9f'{r)f'{s) + -(1 — cos^ 9)f'{r) sin2/(s) 

9 L s 

+ -(1 — cos^ 9)f'{s) sin2/(r) + — (1 + cos^ 9) sin2f(r) sin2/(s) 
r rs 

g 

H-cos 9 sin^ /(r) sin^ /(s) , 


(26) 


(27) 
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where cos0 = x ■ ij, r = \x\ and s = |y|. From (|2^ ) we see that Xf(x)Xf{y) is a function of only three variables, 
Xf{x)Xf{y) = ^ where are immediately found from (p^). They are given explicitly in the 

Appendix for easy reference. We also give in the Appendix the coefficient functions Hn{r, s) which are defined by, 


\ J J / Ti -U 

After carrying out the angular integrations we get an action as a functional of the profile function, 

3 


r = --^ fdrds{rs)'^'^In{r,s)Hn{r,s), 

^ n=0 


(28) 


(29) 


where. 


In{r,s) 



d{cos 9) cos" 9 
(^2 _j_ g 2 _ 2rs cos 9) ’ 


(30) 


are given in the Appendix. 

We were unable to minimize the action with respect to the profile function by analytical methods. We have therefore 
employed a variational method to determine the best profile function approximately. The two important variational 
parameters are the ones that govern the asymptotic behavior of /(r), 


r^oo 


(31) 


The dependence on the parameter a which determines the asymptotics at large distance turns out to be simple. 
Using two trial functions. 


/iW =7^ a , a 
+ x“_ 

/ 2 (r) = 2arctan(p/r 


(32) 

(33) 


and performing the double numerical integration we have found that the action monotonically increases with a. This 
means that the optimal value for this parameter is a = 2, since this is the lowest possible value at which the surface 
terms are non-infinite. 

Having fixed a, we next studied the dependence on the parameter (3 which determines the behavior close to the 
instanton center. The trial functions we used for this purpose are. 


/i(r) = TT 


2//3 


-I- 

/2 (r) = 2 arctan 


{r/pf + {tIpY ) ■ 


(34) 

(35) 


We find numerically the optimal value for the parameter /3 to be /3 = 1.1 for /i and /3 = 1.0 for f^. 
The corresponding values for the action are 


r=1.75-^, (36) 

9 

for /i, and 

r = 1.73^, (37) 

9 

for / 2 . 

We have calculated the first correction AT to the action in the expansion of 5'“^ around Sab for the two profiles. 
We find AF = 0.21-^ for fi and AF = 0.34-^ for / 2 . The correction to the action for both profiles is of order of 
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10%. Note however, that the Ansatz /i is more stable, and indeed after the correction is taken into account has lower 
action than / 2 . 

The corrected (to this order) values for the action are 

o 2 

r = 1.96^, (38) 

9 

for /i, and 

r = 2.07^, (39) 

9 


for / 2 . 

Note that although our Ansatze for the prohle function depend on the instanton size p, the value of the action does 
not depend on it. This is the direct consequence of the dilatational invariance of the effective u-model action eq.(pT|). 
The introduction of the cutoff A strictly speaking breaks the dilatational invariance. This breaking however is very 
small and the invariance is restored in the limit A —> oo. 

We now want to comment on the numerical value of the transition probability obtained in our variational approach. 
Our result eq.(|^) should be compared with the value Sir'^fg'^ for the classical action of a Yang-Mills instanton. So 
the action of an instanton in the variational approach is about two times larger than the standard path integral result 
and the transition probability therefore appears to be much smaller. As explained in Section II it is in fact natural to 
expect this sort of behavior in the Gaussian approximation since the tails of the Gaussians fall off much faster away 
from the minimum than the actual wave function and the overlap is therefore smaller. This is the basic reason for 
the discrepancy in the numerical values of the action between the variational instanton of this section and the weak 
coupling (WKB) instanton of the standard path integral approach. 

It is significant that both the asymptotic behavior and the value of the action for both our variational Ansatze 
for the profile function is very similar to the corresponding results for the Atiyah-Manton Ansatz H- As discussed 
in Section II, the natural identification between the QGD instanton and the instanton in the effective cr-model is 
furnished by the holonomy along all time-lines. 


U{x) = T exp 




For the QCD instanton eq.(|4^) results in the profile function, 


fAM{r) = TT 


I 1 + ^ 


(40) 


(41) 


The asymptotics of this configuration at large distances is f(r) oc l/r^, and near the instanton core /(r) — tt oc t. 
This gives a = 2, /3 = 1 for the Atiyah-Manton configuration, identical with /2 and practically indistinguishable from 
fi. In fact the profile /am is very similar to /i not only asymptotically but also in the whole range 0 < r < oo. In 
Fig. 0 we plot the three profile functions for the same value of the instanton size. It is clear that /i and /am are very 
similar while /2 differs from them somewhat inside the instanton core. We also note that /2 has a slower approach 
to its value at r = 0, which explains why the value of the action for this profile is less stable against the 5“^ — 
corrections. 

We have calculated the cr-model action for the Atiyah-Manton profile numerically and found (including the first 
correction in 3°^^ — expansion) F = 1.97^^. This again is very close to the value we obtain with our variational 
Ansatz eg. (|3^) . This means that even though the value of the transition probability is underestimated in the Gaussian 
approximation, the actual field configurations into which the tunneling is most probable are identified correctly-they 
are precisely the same as in the WKB calculation. 

It should be noted that the exact value of the transition probability due to small size instantons (very weak coupling) 
is not so important, since they do not significantly affect the structure of the vacuum in any case. It would therefore 
be incorrect to declare the Gaussian approximation as a failure from the instanton point of view. One expects the 
Gaussian approximation to do much better at intermediate couplings, since the overlap region contributes significantly 
to the energy and therefore plays important role in the minimization procedure. In the Yang-Mills framework this 
means that there are two types of instantons configurations that we expect to contribute significantly and therefore to 
be better described by the Gaussian approximation. First, these are large instantons for which the coupling constant 
is not small and the action not large. Those are presicely the configurations that cause problems in dilute instanton 
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FIG. 1. The hedgehog profile functions /i, /2 and /am for p = 1. 


gas approximation. The second class is the multi-instanton configurations, where the separation between instantons 
is not too large and the interaction between them is not negligible. These are the configurations of the instanton 
liquid type. As discussed in the previous section, these configurations are in fact responsible for the generation of the 
dynamical scale in the best variational state through the condensation of instantons. The direct discussion of these 
multi-instanton configuration is however beyond the scope in the present paper. The large size instantons, or rather 
the instantons of the size comparable to the dynamical scale 1/M are the subject of the next section. 


IV. MASS CORRECTION AND INSTANTONS OF STABLE SIZE 

We now want to explore the effect of the dynamical mass on the properties of the instantons. Our expectation 
is that the presence of the mass scale stabilizes the size of the instantons and suppresses the instantons of the sizes 
larger than 1/M. 

The simple qualitative argument to this effect is the following. Consider the effective tr-model action for very 
large size instantons. In such a configuration only field modes with small momentum k < M are present. For these 
momenta the action simplifies, and as discussed in the action eq. © becomes the standard local cr-model where 
M plays the role of the ultraviolet cutoff. 

r = J d^U'<{x)d^U{x) (42) 

If the large size instantons are stable at all they should also be present as stable solutions in this local action eq.(p^. 
However this is not the case as can be easily seen by the standard Derrick type scaling argument. Take an arbitrary 
configuration u{x) in the instanton sector and scale all the coordinates by a common factor A. Then obviously 

r[M(Aa;)] = A“^r[it(x)] (43) 

The dependence of the action on A is monotonic and is minimized at X oo. This means that the instantons in the 
local CT-model shrink to the ultraviolet cutoff 1/M. For instantons smaller than the inverse cutoff we can not use the 
local action anymore. However the behavior of these small size instantons is already familiar. We know that when 
the running of the coupling is taken into account, these instantons are pushed to the large size. This is the infrared 
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problem of large instantons we alluded to earlier. In our variational state the coupling constant stops running at the 
scale M. The picture is therefore very simple. The small size instantons are pushed to larger size by the effect of the 
coupling constant, while the large size instantons are pushed to smaller size by the effect of the local tr-model scaling. 
It is therefore clear that the size will be stabilized somewhere in the vicinity of p ^ I/M. 

To confirm this picture we now turn to the numerical evaluation of the instanton action. We write the inverse 
propagator as, 

G-i(fc) =G-i(fc) + AG-i(fc), (44) 


where 


Go\k) = \k\ (45) 

AG-\k) = 9{M - |fc|) (M - |A:|) 

Since the action eq.([ll|) is linear in G“^, the Gq^ term will lead to the classical action Tq calculated in the previous 
section. As discussed above, Tq is scale invariant and therefore does not depend on the size of the instanton. In 
coordinate space, the second term in eq.(^ reads, 

AG"^(|a; - y\) = - ^|4 (cosM\x - y\ - 1 + ^K mM\x - y\^ . (46) 

The mass correction term in the propagator gives rise to the correction in the action 

AT = i I d^xd^yXnx)AG-\x - y)Xt{y). (47) 


We now compute the dependence of AF on the size of the instanton for the three profile functions considered above, 
/i, /2 and /am- As before, the angular integrations can be performed exactly. We obtain. 


AF = -- 


where are as before and defined by 


drds{rsY E kn {t , (r, s), 


(48) 


n—0 


In = 


^ (i(cos 0) cos" 0 
_i (r^ + _ 2rscos0)^ 


cosMja; — y| — 1 + ^ sinMja; — y\ \ , 


(49) 


are given in the Appendix. 

We calculated numerically the dependence of the classical action^ on the size of the instanton for the three profiles 
(Q), (^^ and (|^). We find that the action increases monotonically with the instanton size. The dependence of the 
classical action on the size for the three anzatse is shown in Fig. ||. As discussed in the beginning of this section the 
classical action favors small size instantons in all three Ansatze. 

We now must include the effect of the running coupling constant. This is equivalent to calculating the one loop 
correction around the instanton background. Rather than performing this technically nontrivial calculation we will 
use the results of |^]. We will limit ourselves to the SU{2) theory in the following. The /3-function in the nonlinear 
CT-model was found in as 


= -^ 8 - ( 50 ) 

This is slightly different from the complete one loop Yang-Mills /3-function. The origin of this discrepancy was studied 
in and is well understood now. It is due to the fact that the present variational calculation omits the screening 
contributions of the transverse gluons. As discussed in [Q the variational state can be modihed to yield an exact 


^In this section we limit ourselves to the approximation 5'“*’ = 5“*’. The corrections to this approximation as we saw in the 
previous section are fairly small and we will not explore them here. 
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FIG. 2. Dependence of the classical action on the size of the instanton for the profiles /i, /2 and /am- 
in units of 


The action is given 


/3-function. This point is not essential for our analysis and we will not dwell on it any further except noting that the 
use of either eq.(|^) or the exact one loop expression in the following leads to the same results. 

Eq.(pO|) is valid at distances smaller than 1/M. The coupling constant at these distances therefore scales according 
to 


87r2 
> M) 


= 81n 



(51) 


At distances larger than M the running of the coupling constant stops and it tends to a constant value, 

87r2 _ o 1 ^ 

g^{g<M) \J^) ■ 


(52) 


with a, a numerical constant of order one. The precise interpolation between the two regime is unimportant and we 
will use the following simple interpolating expression 


87r^ 

g‘^{g) 


= 4 In 


M2 

^'^QCD 



(53) 


Since the exact value of the constant a is not known we will present our results for several values of order one. When 
evaluating the action of the instanton of the size p we obviously must take p, = 1/p. 

The running of the coupling constant is not the only logarithmic effect at one loop order. In addition one has to 
take into account the path integral measure over the three translational zero modes and the size of the instanton, 

dxidx2dx3dpp~'^, (54) 

which contributes an extra term to the action, 

Tmeasure — 4 log p. (^^) 


Collecting all the one-loop logarithmic contributions together we calculate the corrected action as a function of the 
size of the instanton as 
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FIG. 3. The classical action plus one-loop logarithmic corrections for the profile fi as a function of the instanton size. The 
different curves correspond to different values of the parameter a = 0,..., 2, and the bullets mark the minimum of each curve. 


a 


fi 



/2 



/am 


Pmin 

r(pmin) 

r"(Pmin) 

Pmin 

r(pmin) 

r"(Pmin) 

Pmin 

r(pmin) 

r"(pn,in) 

1/3 

1.5 

32 

6.9 

0.9 

36 

23 

1.6 

32 

6.3 

1/2 

1.4 

33 

7.9 

0.8 

37 

26 

1.5 

33 

7.0 

1 

1.1 

36 

9.7 

0.7 

39 

33 

1.2 

36 

8.1 

2 

0.9 

39 

12.0 

0.6 

41 

41 

1.0 

39 

9.7 


TABLE I. Stable size pmin for the soliton in units of 1/M = 0.75GeV the value of the minimum action r(pmin) and the 


curvature at the minimum r"(pniin), for the three profiles, for different values of the parameter a. 


r(5^(i/p)) + rmeasure- (56) 

The result for the three different proHles we have considered are shown in Figs. and ||, for different values of 
the parameter a. 

In all cases the action has a minimum at a size of the order of 1/M. The values for the size of the instanton, its 
action and the curvature at the minimum are given in Table || for different values of a. The results for the Atiyah- 
Manton Ansatz and our function fi are practically indistinguishable. Varying the parameter a between 1/2 and 2 
the size of the stable instanton varies between 1.5/M and 1/M. The Ansatz /2 leads to the instanton size smaller by 
a factor of 1.5 — 2. We consider however the former estimate to be better, since the Ansatz gives a considerably 
larger value of the action (see Fig. 2) and is therefore not a very good choice for the instanton “valley” conhguration. 

The numerical analysis of this section confirms therefore our expectations. The presence of the dynamical mass sta¬ 
bilizes the size of the instantons at the value p = ■ The large instanton problem therefore finds a nonperturbative 

solution in the framework of the Gaussian variational approximation. 


V. DISCUSSION 

In this paper we have discussed how the instanton transitions appear in the framework of the gauge invariant 
Gaussian approximation of |^] . The relative gauge rotation between two gauge field configurations between which the 
tunneling transition is most probable is determined from the saddle point equation of the effective nonlinear a-model. 
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FIG. 4. The classical action plus one-loop logarithmic correction for the profile /2 as a function of the instanton size. The 
different curves correspond to different values of the parameter a = 0,..., 2, and the bullets mark the minimum of each curve. 
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FIG. 5. The classical action plus one-loop logarithmic correction for the profile /am as a function of the instanton size. The 
different curves correspond to different values of the parameter a = 0,..., 2, and the bullets mark the minimum of each curve. 
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We found that the relative gauge transformations for most probable tunneling transitions are very similar to the ones 
found in the standard path integral instanton approach. 

The value of the logarithm of the tunneling probability for small size instantons however turns out to be larger by 
about a factor of two in the variational vacuum. This is understandable since the tail of a Gaussian wave function 
decreases faster and therefore the overlap integral of two such wave functions is smaller than of the semi-classical wave 
functions. 

Our main result concerns the effect of the dynamical mass scale that characterizes the variational vacuum. We 
find that the presence of this scale stabilizes the size of the instanton. The instantons of the size larger than the 
inverse of this scale are suppressed. When account is taken of the proper running of the strong coupling constant, 
the integral over the instanton sizes has a saddle point. This saddle point determines the most likely size of the 
instanton as p = 1 — 1.5/M. The large size instanton infrared problem which plagues the standard dilute instanton 
gas approximation is thereby removed in the variational vacuum due to the presence of the dynamical nonperturbative 
scale. 

An interesting point is that the instanton action for the most likely size instanton is pretty large, the numerical 
value being around 35 (see Table This means that the configurations with small number of instantons and anti- 
instantons are not important energetically. Nevertheless, the generation of the dynamical mass itself within the 
framework of our variational approximation is due to the instanton condensation. This is so since the effective a- 
model is in the disordered phase at the value of the best variational parameter M. This suggests that the important 
type of configurations are the ones that contain many instantons and anti-instantons. The small fugacity factor in 
these configurations can be overcome by a large entropy and also by the interaction between instantons and anti- 
instantons. This interaction is known to be attractive for some relative color orientation. The interaction also is quite 
long range, since the instanton profile function away from the instanton center decreases only as the second power of 
the distance. It is therefore entirely possible that the important configurations for the energy minimization are of the 
type of the instanton liquid-namely those having large number of instantons which although fairly dilute nevertheless 
feel each others presence strongly due to the long range interaction. It is in fact interesting to note that the most 
likely instanton size we found here is consistent with the average size of the instantons in the instanton liquid model of 
|^,Q. For the case of SU{2), the average instanton size, in units of the gluon condensate obtained in the instanton 
liquid model, turns out to be 

P ~ .4. 

In our case, taking the value of the gluon condensate obtained in the variational approach |^, we find. 

The relation of the variational approach with the instanton liquid model is a very interesting open question and 
warrants further study. 
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APPENDIX: 


The finite action for the case of the mass scale M = 0 is given in ( p6[ ) . To obtain the explicit action for the profile 
of the hedgehog field (^) it is necessary to calculate the Laplacian over v = x — y oi the product of the right currents 
A“ at the points x and y, which has the form 


= 4^ cos” s). 

^ n=0 


The expression for Hn are obtained from (p^. 


(Al) 
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Ho = — sin2/(r) sin2/(s) + - sin2/(r)/'(s) + - sin2/(s)/'(r), 


(A2) 


rs 


r 



(A3) 


H 2 = — (sin2/(r) - 2r/'(r)) (sin2/(s) - 2s/'(s)) . 
rs 


(A4) 


After calculating the Laplacian over v we obtain 



(A5) 


From Eq. (|^) the expressions for Hnir, s) can be found: 


2r^s^Hoir, s) = 8rs sin^ /(r) sin^ /(s) + sin2 /(r) sin2 /(s) 

— 16 s cos/(r) sin/(r) sin^ f(s) f'{r) — 2 sin2 /(s) f(r) 

—2 sin 2 /(r) sin 2 /(s) f {r)^ — 2r^s sin 2 f{r) f'{s) 

+4r^ s f{r) f{s) - 4r^ sin2 /(s) /'(r) f'{sf + 2r^ sin2 /(s) fir) 

+r^ cos2/(r) sin2/(s) /"(r)+2r^s^ cos 2/(s)/'(r) /"(s) 

+r^ sin2/(s)+ (r ^ s), (A6) 

2r^s^Hi{r, s) = —8r^ sin^ /(r) sin^ /(s) + 3r s sin2 /(r) sin2 /(s) 

—4r^ s sin2 /(s) /'(r) — 2r^ s cos2 /(r) sin2 /(s) /^(r) 

+8 cos /(r)^ sin^ /(s) f'{r)^ + 8 cos 2 /(r) /'(r) f'{s) 

—4r^ cos2 /(r) cos2 /(s) /'(r) /'(s) — 8 sin^ /(r) sin^ /(s) /'{sf 

+8 cos/(r) sin/(r) sin^ /(s) f''{r) + 6 s sin2 /(s) /”(r) 

-8r^s^ fis)f"ir) - 4r^s^ cos2/(s)/'(s)/"(r) + (r ^ s), (A7) 

2r^s^H2{r, s) = —32 r s sin^ /(r) sin^ /(s) — 3 sin2 /(r) sin 2 /(s) 

+32 s cos /(r) sin /(r) sin^ /(s) /'(r) + 6 sin 2 f{s) f'{r) 

—2 sin2 /(r) sin2 /(s) f {r)^ + 6r^ s sin2 /(r) /'(s) 

—12s/'(r)/'(s) — 32cos/(r) cos/(s) sin/(r) sin/(s)/'(r)/'(s) 

+4r^ sin2 /(r) f{r)^ f'{s) — 2r^s^ sin2 /(s) f"{r) 

+r^ cos2/(r) sin2/(s) /"(r)+4r^ f'{s) f''{r) 

—2r^ cos2 /(s) /^r) /"(s) — sin2 /(s) 

+2r^ /'(s)+ (r s), (A8) 

2r^s^Ho{f, s) = —9rs sin2/(r) sin2 /(s) + 12s sin2 /(s) /^r) 

+6r^ s cos2 /(r) sin2 /(s) /'(r) — 16 /^(’') /^(s) 

—8 cos 2 /(r) /'(r) /'(s) — 4 cos 2 /(r) cos 2 /(s) /'(r) /'(s) 

—6 s sin 2 /(s) f'{r) +8r^ f'{s) f''{r) + 4 cos 2 /(s) /'(s) 

—4r^ /"(r)/"(s) + (r s). (A9) 


The expressions for the angular integrals (^) appearing in the action (p^ are the following: 



(AlO) 
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The corresponding integrals (^9|) for the case of M 7 ^ 0 read: 

1 / sin^(t/ 2 ) Y 


/n = 


h = 


2rs 


J 

'-I0 + 


M\r — s\ 

1 / cost 


f ^ f 
^2 — —— 


2j^2g2 

1 


X M(r+s) 

- ci(t) + ln(t/2) I 

/ M\r — s\ 


[(—4 + — s^) cost 


+ 2 (r 2 +s 2 )(ci(t) - In (t/ 2 )) + t(t - 4sint)]^[^;^"l , 
where ci(t) is the Cosine-integral function 

cost 




-dt. 


(All) 

(A12) 

(A13) 


(A14) 
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